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Abstract: The IceCube neutrino spectrum shows a flux which falls of as E ~‘^ for sub 
PeV energies but there are no neutrino events observed above ~ 3 PeV. In particular the 
Glashow resonance expected at 6.3 PeV is not seen. We examine a Planck scale Lorentz 
violation as a mechanism for explaining the cutoff of observed neutrino energies around a 
few PeV. By choosing the one free parameter the cutoff in neutrino energy can be chosen 
to be between 2 and 6.3 PeV. We assume that neutrinos (antineutrinos) have a dispersion 
relation E"^ = — {^ 3 /Mpi) p^, and And that both 7 r+ and 7 r“ decays are suppressed 

at neutrino energies of order of few PeV. We find that the /r“ decay being a two-neutrino 
process is enhanced, whereas decay is suppressed. The K~^ —>■ is also suppressed 

with a cutoff neutrino energy of same order of magnitude, whereas K~ —)■ 7r®e“Pe is 
enhanced. The n —>■ p~^e~i'f, decay is suppressed (while the h ^ p~e'^i'e is enhanced). 
This means that the Ug expected from n decay arising from p + 7 —A — 7 r'’'n reaction 
will not be seen. This can explain the lack of Glashow resonance events at IceCube. If no 
Glashow resonance events are seen in the future then the Lorentz violation can be a viable 
explanation for the IceCube observations at PeV energies. 
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1 Introduction 

IceCube collaboration has observed the neutrinos of very high energy going to beyond 2.6 
PeV [1-4]. The IceCube data in the energy range 60 TeV to ~ 3 PeV is consistent with 
E~‘^ neutrino spectrum following E^dNy/dEy ~ 1.2 x 10“^ GeVcm“^s“^sr“^ [2, 3]. A 
neutrino spectrum sharper than E~^'^ does not give a good fit to the data [3]. There are 
no neutrino events observed above ~ 3 PeV. 

In particular, there is no evidence of Glashow resonance [5], + e~ —)• W~ —)• shower, 
which is expected at E = 6.3 PeV. Glashow resonance gives rise to an enhanced cross- 
section for Vf, at resonance energy E = M^/2me = 6.3 PeV, which increases the detection 
rate of Ug + 9ehy a factor of ~ 10 [2]. This implies that at least three events should have 
been observed at Glashow resonance, but none were. 

The Glashow resonance gives rise to multiple energy peaks at different energies [6]. The 
first one is at 6.3 PeV and others lie at the Eyis = E — Ex, where Ex is the energy in 
the W decay, which does not contribute to the visible shower [7]. The decay of W into 
hadrons goes as VP —)• qq, giving rise to a peak at 6.3 PeV, while decay into leptons goes as 
W Dl, which means W boson will lose half of its energy and so a second peak at 3.2 PeV 
is expected. In case of r lepton in the final state, a further decay takes place producing a 
neutrino and thus a third peak at 1.6 PeV. The events observed by IceGube [1-4] between 1 
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PeV to ~ 3 PeV range may be associated with the second (leptonic decay of W) and third 
peak (r decay), but non-appearance of Glashow resonance hadronic shower from W ^ qq 
at 6.3 PeV (dominant peak) makes this idea less attractive. The non observation of the 
expected signature of Glashow resonance in IceCube data indicates a cutoff of neutrino 
energies between 2-6.3 PeV [7, 8]. 

In this paper, we propose a mechanism which can explain why neutrinos above a certain 
energy may be suppressed in the astrophysical production processes like tt —)• K —)• 

etc. We assume that Lorentz violating higher dimensional operators [9, 10] give rise 
to a modified dispersion relation for the neutrinos (antineutrinos) of the form E'^ = p'^ + 
m‘1 — p^ with n > 2. Depending on the sign of ^n, the neutrinos (antineutrinos) 

can be either superluminal {^n < 0) or subluminal > 0). For the superluminal case, 
it has been shown [11, 12] that the presence of the extra terms in the dispersion results 
in a suppression of tt and K decay widths. The phase space suppression for both the 
subluminal and superluminal dispersions for meson decay and the Gerenkov process v —^ 
iye~^e~ has been noticed in [9, 13-16] with limits on Lorentz violation parameters from 
IceGube events. A comprehensive listing of Lorentz and CPT violating operators and their 
experimental constraints is given in [17]. In this paper, we calculate the vr, K, p and n decay 
processes in a fixed frame (the frame chosen being the one in which the GMBR is isotropic; 
although the Earth moves at a speed vEarth ~ 300 km/sec with respect to the GMBR, the 
Lorentz correction to the neutrino energy is small as f^Earth ~ 10~^)) where the neutrinos 
(antineutrinos) dispersion relation is — {^^/Mpi) p^ [10, 18-20]. We will have 

^3 > 0 for neutrinos and ^3 < 0 for antineutrinos. In the 7 r+ decay, we hnd that the |Mp is 
suppressed at neutrino energy E^,, where — m?^ ~ {^s/Mpi) p^. This implies that for the 
leading order Planck suppression (n = 3) taking ^3 ~ 0.05, the decay is suppressed at 
Ey ~ 1.3 PeV. Similarly decay will be cutoff at Ey ^ 2 PeV with rn\—rn?^ ~ {^ 2 ,/^pi)p^ 
and neutron decay will be cutoff for p, where [mn — nip)'^ ~ (? 3 /Afpz)p^, which is lower than 
the Glashow resonance energy. For the tt~ decay the \M\‘^ is enhanced but the phase space 
is suppressed and therefor vr” —)• is also suppressed. In the case of p~ e~9e^p 

decay, |Mp is enhanced whereas the phase space suppression is not significant, so the p~ 
decay rate is enhanced (while —>■ decay rate is suppressed). This enhancement 

is significant at p~ energies ~ 2 PeV but since the primary source of p~ is vr” decay which 
is already cutoff, there will be no observable effect of this enhancement in the neutrino 
spectrum seen at IceGube. Neutrinos from K~ —)• p~9p_ and —)• decays will 

be cutoff at slightly higher energies. Radiative vr^ decay with a single neutrino in the 
outgoing state are also suppressed. The three body kaon decay rate are determined by the ^3 
dependence of |Mp and we find that K~^ —^ p'^Vp decay is suppressed but K~ —)• -k^ p~Vp 

decay is enhanced. Neutron beta decay n —>■ p'''e”7'e gets suppressed in the same way as p'^ 
decay. If the source of T'e is neutron beta-decay [21] then the mechanism proposed in this 
paper can be used for explaining the absence of Glashow resonance [5] at IceGube. The 
value of i^s/Mpi) rsj 0.05 used in this paper to explain the cutoff in PeV neutrinos is 

much smaller than the bound on the dimension-five coefficient, (aif^)oo < 3.5 x 10”^*^ GeV”^ 
from SN1987A dispersion [13]. 
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The rest of the paper is organized as follows. In section 2, we calculate the leptonic decay 
widths of pions and kaons using modified dispersion relation of neutrino and compare them 
with their standard model counterparts. In section 3 we study —)■ K~^ —>• 7r^e'’'t'e 

and n —)• p'^e~i>e processes with modified neutrino dispersion. We give our conclusion in 
section 4. 


2 Two body decays 

2.1 Neutrino velocity with modified dispersion 


To calculate the decay widths of pion, kaon and muon, we use the following dispersion 
relation, 

t7i2 _ 2 I 2 yi / n -1 \ 

E =P (2.1) 

which is motivated by Lorentz violating higher dimensional operators [9, 10]. We will take 
> 0 for neutrinos and < 0 for antineutrinos. We use this modified dispersion relation 
to get the neutrino (antineutrino) velocity, which becomes 


- ^ - 1 _ ^ ~ ^ 'Cn 
dp 2 


( 2 . 2 ) 


This is clear from eq.(2.2) that we have a subluminal neutrinos and superluminal an¬ 
tineutrinos. In this paper, we will consider the leading order Planck suppressed dispersion 
relation E‘^ = p'^+ —{^ 3 /Mpi) p^ to compute the primary decay processes which produce 
neutrinos and antineutrinos. In appendix. A, we obtained modified dispersion relations for 
neutrinos and antineutrinos using dimension 5 operator. 


2.2 7 r+ p+iyp 

We calculate the pion decay width using the modified dispersion relation for neutrino by 
taking n = 3 case. The amplitude calculation of pion decay process TT~^{q) —>• pi'^{p)vp{k) 
gives, 

M = fy^Vud g^u(A:)^7^(l - 75 )u(p) (2.3) 

where /^r = /(m^) is a constant factor and Vud is the CKM matrix element. The spin 
averaged amplitude squared is. 


|M|2 


2Glfl\V^d?mlF{k) 


2 2 
m^-rrip 





(2.4) 


where ^3 = Cs/Mpi and the F{k) factor comes from the modified spinor relation of neutrino, 
as described in eq.(B.9). The decay width of pion is then given by, 


r = 


Gy^Kdl^mlFik) 




X 




k‘^ dk d cos 9 

Eu\J\q- + rrif^ 

2 


-AEu. -Fyy + J\q- A :|2 ml) 


m. 


+ 2 


m 




(2.5) 
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Figure 1. The ratio F/Tsm for 7 r+ —>■ and 7 r“ —^ processes in Lorentz invariance 

violating framework to its standard model prediction for superluminal (^3 < 0 ) and subluminal 
’-'fj. (^3 > 0) final states as a function of pion momentum We considered ^3 = ±1.3 x 10“^ for 
corresponding processes. 


after using = F{k)k, and writing |pl = \q — fcp = k'^ — 2kqcos6, our expression of 

eq.(2.5) takes the following form 


r = 


f kdkdcosO 


/ 


SttEtt 


q — k\^ ± m? 


mz 


- E^ + J\q- A:|2 ± m^) 


mt-ml-esk^(^ + 2 


from the argument of the delta function in eq.(2.6), we have 


which gives, 


cos 9 = 


^J\q-k\^+ml = E^-E,^ 

(^2 _ ^2 + 2 E^k - E^k% + k%) 


2 kq 


We reduce the 6 function in E^ to a d function in cos 6 by taking, 


d 


d cos 9 


{Ey. -E^ + J\q- A:|2 ± ml) 


kq 


^kP‘ + q^ — 2kq cos 9 ± 


( 2 . 6 ) 


(2.7) 


( 2 . 8 ) 


(2.9) 
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Figure 2. The ratio T/hsM of tt’’' —> process in Lorentz invariance violating framework to 

its standard model prediction for subluminal neutrino (^3 > 0 ) as a function of neutrino energy 
kmax with different values of ^ 3 . 


and substituting in eq.(2.6). We get the pion decay width, 


r = 


SttEt, 


f dk 

J q 


^ + 2 


mf 


( 2 . 10 ) 


We solve the integration in the limits of k, which are hxed by taking cos0 = ±1 in eq.(2.8), 


k 


max 


+ ^3^maxi^Tr kmax) 

2{E^ - q) 


( 2 . 11 ) 


h ■ — 


m: 


^'^kmini^T^ kmin) 


2(Etj + q) 


( 2 . 12 ) 


solving these equations numerically, we get the allowed limits of neutrino momentum. We 
solve eq.( 2 . 10 ) and then compare our result with the standard model result of pion decay 
in a moving frame, which is 


r5Af(vr fJ.1^) 


SttE-jt 



(2.13) 


We compute the pion decay rate numerically for superluminal Ue (^3 < 0) and subluminal 
L'e (Cs > 0 ) final states and obtain the following : 
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Figure 3. The maximum neutrino energy, k^ax as a function of Lorentz invariance violation 
parameter ^3. 


• For subluminal neutrino final state (.^3 > 0), the allowed phase space (eq.2.11-eq.2.12) 
goes up but the |Mp (eq.2.4) is suppressed. There is a net suppression in r(7r''’ 

as shown in Fig.(l) for ^3 = 1.3 x 10“^. 

• For superluminal antineutrino final state (.^3 < 0), the phase space (eq.2.11-eq.2.12) 

is suppressed but the |Mp is enhanced. The net effect however is a suppression in 
the F(7r“ —>■ for this case also [11], as shown in Fig.(l) for ^3 = —1.3 x 10“^. 

In Fig. (2), for the process 7r+ —)• we show the maximum neutrino energy for different 

values of ^3 using the solution for q in terms of kmax and kmin from eq.(2.11-2.12) in 
eq.(2.10). We see that for ^3 = 5.0 x 10“^, the neutrino spectrum cutoff at kmax = 1-3 
PeV. The upper limit of observed neutrino energy provides bound on the Lorentz invariance 
violation parameter ^3. In Fig. (3), we show the maximum neutrino energy kmax, as a 
function of Lorentz invariance violation parameter ^3. This is clear from Fig. (3) that kmax 
goes down as ^3 increases. 

2.3 K+ 

In the similar way like pion decay, we calculate the kaon decay width for the process 
using the modified dispersion relation for neutrinos by taking n = 3 
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Figure 4. The ratio T/Tsm of K'^ —?► process in Lorentz invariance violating framework to 

its standard model prediction for subluminal neutrino (^3 > 0 ) as a function of neutrino energy 
kmax with different values of ^ 3 . 


case. We get the kaon decay width, 


r = 


Glfl\V, 


us\^mi 


SttEk 


/ 


q 


2 2 

rriK-m^ 





(2.14) 


In the same way like pion, we solve the integration in the limits of k by taking cos 0 = ±1 
which gives, 


k 


max 


- ml + ^'^kl^^^{EK - kmax) 
2{Ek — q) 


(2.15) 


kmj.n — 




li^K kmin) 


2 {Ek + q) 


(2.16) 


solving these equations numerically, we get the allowed limits of neutrino momentum. We 
solve eq.(2.14) and then compare our result with the standard model result of kaon decay 
in a moving frame, which is 


^SmIK nu) 


G],fK\Vus\‘^rnlm 

8 ttEk 


2 

K 



(2.17) 


In Fig. (4), we show the maximum neutrino energy for different values of ^3 using the 
solution for q in terms of kmax and kmin from eq. (2.15-2.16) in eq.(2.14). We see that for 
^3 = 5.0 X 10“^ the neutrino spectrum cutoff at kmax = 2 PeV. 
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3 Three body decays 

3.1 e~Vei'fj_ 

We compute the muon decay width with subluminal neutrino and superluminal anti¬ 
neutrino in the final state, assuming the dispersion relation for the neutrino (antineutrino), 
E‘1 = k‘^ — where ^3 > 0 and .^3 < 0 correspond to subluminal neutrino and superlu¬ 
minal antineutrino respectively. We assume identical ^3 for all the species of v (and z^) to 
avoid an extra source for neutrino oscillations which is not observed [16, 22]. The amplitude 
for the process ^~{p) —)• e~ is given as, 

M = ^[u{k')'y^{l - -/5)v{k)][u{p')-fi,{l - 'IbMp)] (3.1) 

where Gi? is the Fermi constant. After squaring amplitude and solve it using trace tech¬ 
nology, we get the spin averaged amplitude, 


\M\^ = 64Gj,{p ■ k){p' ■ k') 


(3.2) 


The decay width of muon is, 


dr= - ^-^{27r)'^5'^{p - p' - k' - k) 

(27r)32E,^ (27r)32Se (2vr)32Sp^ 2E/ ^ ^ ^ 

using the squared amplitude from eq.(3.2), we get 
32 d^k' d^p' d?k ^4 


8(27r)3.F^ Ee E^^ Ep^ 
First we write eq.(3.4) as, 

on r 

rp 


6 ^{p — p — k — k){p ■ k){p ■ k ) 


F = 


32 Gl r d^k 


8 { 2 TT)^Ef, 


j '^p'^k'^Ia,fi{p - k') 


where 


^I3{p-k') = J 


d^k d^p' 
Er. E,,„ 


S^{p-p' -k' - k)kap'i3 


(3.3) 

(3.4) 

(3.5) 

(3.6) 


and then to find out Iai 3 {p — k'), we use the generic phase space integral formula, 
d^p 


^ ol (5 




,2 ^ ^ / 2 ^ ^ J ^{k-p- q)paqp = -^{k^[k^ - (jni-m2f] 

y/m^+p-p^/mf + q-q ^ 2 k 


[k^ - (mi m 2 ) ] 5 a /3 + 2[A: -h k^{mi + m^) - 2{m\ - m^) ]kakp) 


where 


27r 


I = ■p-\/[A:2 - (mi - m2)^][A;2 - (mi -h m2Y]. 


(3.7) 


(3.8) 


Applying this to our scenario by putting m\ = m?^ = "^2 = taking 

k = p'/2 ~ p/4, we find 

7 CsP^ 


vr 


Iap{p-k') = - 


1 + 


64 {p - k'Y 

u /^2 , ^ 


(3.9) 


[{p — k) + -^^sP ]9al3 + 2 




1 - 


7 ggP' 

Mip-k'y 


{p-k')a{p-k')fi 
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Figure 5. The ratio T/hsM for —>■ and pT —>■ processes in Lorentz invariance 

violating framework to its standard model prediction for superluminal antineutrino (^3 < 0 ) and 
subluminal neutrino (^3 > 0) final states as a function of muon momentum . Here we considered 
^3 = ±5.0 X 10-2. 


after contracting lap with the muon and electron momentums which respectively are p and 
k' , we get 


lap {p- k') = 


TT 


where, 




1 ± 


7 

64 {p - k'Y 

7 


(3.10) 


[{p-k'f +—i'^p\p-k')+ 2 


7 

64 {p - k'Y 


1 - 


{p ■ p — p ■ k'){p ■ k' — k' ■ k')^ 


p-p = m^ 

k' ■ k' = ml ^ 0 


p ■ k' = k' {Efj, — p cos 6) 

{p — k')"^ = — 2k' {E^ — pcosO). 


(3.11) 


The decay width from eq.(3.5) can be written as, 

32Gl ( 27 r) , pml/ 2 {E^-kcose) 

r = „ ^ — / a cos 6 


8(27r)® El 


At J-i 




k'dk'p^k'^Iap 


(3.12) 
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after solving it, we finally get, 


r = 


l927r^Ef, 





(3.13) 


We compare our result with the standard model prediction of muon decay in a moving 
frame, which is 

rsM(^ ^ (3.14) 

We compute the muon decay rate for subluminal neutrino (^3 > 0) and superluminal 
antineutrino (^3 < 0 ) and obtain the following: 


• The decay rate of the process r(/i —>■ e is enhanced, as shown in Fig. (5) for 

^3 = ±5.0 X 10-2. 

• The decay rate of the process r(/x+ ^ e^VeVfj) is reduced, as shown in Fig. (5) for 
^3 = ±5.0 X 10-2. 


3.2 K+ 7r°e+i^e 

We also calculate 3-body kaon decay width using the modified dispersion relation for neu¬ 
trino by taking n = 3 case. The amplitude calculation of kaon decay process K^{pk) —)• 
7r^{pT,)e+{pe)i^e{Pu) gives. 


|M|2 = 16Gp\Vus\^fl[m^K{pK-Pu+P7z-Pu)-‘i{pK-Pu){PK-PTT)-‘2{pK-Pu){PK-Pu)-mK^'^pl] 

(3.15) 

where /+ is the kaon form factor. The Decay width of kaon is, 


dr = 


which gives, 


d^p-j, (fp^^ (fpe |M|2 

(27r)32F;^ (27r)32F;^^ (27r)32F;e 2Ek 


(27r)^(5^(px -P-K- Pu, - Pe) 


F ~ 


Gp\Vus\ flrn\ 


m\ ( 1 - 


8 m 2 


m 


K 


- ( 1 - 


mZ 


m 


K/ -i 


(3.16) 


(3.17) 


It is clear from eq.(3.17) that the K~^[K ) decay rate goes down (up) as kaon momentum 
Pk increases, which is shown in Fig. ( 6 ) for ^3 = ±5.0 x IO- 2 . 


3.3 n —>■ p~^e Ve 

In the similar way like muon decay, we also calculate the neutron beta decay width using 
the modified dispersion relation for antineutrino. The spin averaged amplitude squared for 
the neutron decay process n{p) p'^{k)e~ {k')r'e{p') comes, 

\M\^ = 64Gl{p ■ p'){k ■ k') (3.18) 


using eq.(3.18), we get the following decay width of neutron, 


dF = 


32 Gj^ d^kd?k'd?p' 
Ee Ep^ 


8(27r)5£'„ Ep 


5^{p — k — k' — p'){p ■ p'){k ■ k') 


(3.19) 
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Figure 6. The ratio F/Tsm for —>■ 7r°e+z^e and K~ —>■ Tr^e~Pe processes in Lorentz invariance 

violating framework to its standard model prediction for superluminal Pe (^3 < 0 ) and subluminal 
iZe (^3 > 0) final states as a function of kaon momentum pK- We considered ^3 = ±5.0 x 10“^ for 
corresponding processes. 


we solve eq.(3.19) in the similar way like muon decay using generic phase space integral 
formula (eq.3.7). Then we solve the final integral over the electron energy, for which 
the minimum energy is the rest energy nie of the electron while the maximum energy is 
approximately, 


Er, 


rrin - nir. 


(3.20) 


which finally gives. 


G]p{mn - mpf-nir, 
Ibir^En 


{rUn - mpf - 


(3.21) 


For '^3 = 0.05 the neutron decay width goes down at neutrino momentum p ~ 0.1 PeV. 
This implies that antineutrino production from neutron decay will be suppressed and so in 
our model, it is also possible to explain the absence of Glashow resonance [5]. The decay 
rate of the charge conjugate process n —)■ Penh's is enhanced, but since only neutrons are 
produced in the pTy—)-A—)-n± Tr~^ processes at the source, the enhanced decay of n is 
not relevant to the IceCube events. 
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4 Conclusion 


In this paper we provide a mechanism by which one can account for the lack of an¬ 
tineutrino events at Glashow resonance (6.3 PeV) at IceCube. We show that if the 
neutrino (antineurino) dispersion is modified by leading order Planck scale suppression 

= — {^■i/Mpi)p^ (where ^3 > 0 correspond to neutrinos and ^3 < 0 correspond to an¬ 

tineutrino) , then there is a suppression of the 7 r+ decay width and corresponding neutrinos 
will be cutoff at energies Ey = 1.3 PeV (with ^3 = 0.05). The neutrinos from Kaon decay 
K+ —>• will be cutoff at 2 PeV. 

• Three body decays like pL~ e~r'ei^fi and K~ —)• Tr^e~i'e get enhanced due to different 
^3 dependence in their |Mp, whereas three body decay widths of and K+ get 
suppressed. 

• Neutron decay n ^ p'^e~v'e gets suppressed in the similar way as decay. So if the 
source of i>e is neutron beta-decay then the mechanism proposed in this paper can 
be used to explain the absence of Glashow resonance at IceCube. 

• Radiative three body decays like 7 r=*= —)• e^z /7 and —)• //^z /7 are factorized to the 
|Mp for two body decays —)• e^u and —)• times ctem [23, 24] and these are 
also suppressed like two body decay processes. 

The enhancement in p~ decay will be significant at muon energies of 2 PeV and if the 
primary source of p~ is 7 r“ decay then there will be no observable consequence of this in 
IceCube events. However such enhancement of the p~ decay rate would be observable for 
p~ produced not from 7 r“ decay but e.g via pair production e.g in e“*"e“ —)■ p'^p~. The 
precise numerical values depend on the choice of the parameter ^ 3 , but obviously a cutoff 
between ~ 3 PeV and 6.3 PeV can be easily obtained in this model. We conclude that 
if neutrinos at Glashow resonance energies are not observed at IceCube then explanations 
in terms of new physics such as Lorentz violating modihed neutrino dispersion relation 
become attractive. The fact that neutron decay into p + e + i?e is suppressed has the 
following implications. The conventional vr/K decay neutrinos from astrophysical sources 
have cutoff in the range of ~ 3 PeV. However the B-Z neutrinos which arise in GZK process 
have two components [25], the higher energy neutrinos from vr/K will be more suppressed 
compared to the lower energy n decay to Ue- But both components of GZK process will be 
suppressed at Ey > 3 PeV. 
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A Dispersion Relation 


The cubic dispersion relation we used for neutrinos and antineutrinos can be obtained from 
the dimension 5 operator [9, 10], 

^LV = + (A.l) 

Mpi 

where is a fixed four vector that specifies the preferred frame. Both the vector and 
axial-vector terms in eq.(A.l) are CPT violating in addition to being Lorentz violating. 
The Lagrangian gives the equation of motion, 

+ (A. 2 ) 

where we have taken E ^ m. This leads to the following dispersion relation for left and 
right handed particles ' 0 , 

3 

E^=p^ + 2{m±V2)^ (A.3) 

Mpi 

where -I- and — signs correspond to V’J? and V’L respectively. Now taking the charge conju¬ 
gation of eq.(A.l), we find 

^LV = (A.4) 

where we used charge conjugation properties viz. = — 7 ^ and C'“^ 7 ^ 75 C' = 7 ^x 75 . 

The operator (eq.A.4) gives the following dispersion relation for left and right handed 
antiparticle 

3 

_^2(-r/i ± 72 )^ (A.5) 

where the -|- sign is for V’/j and — sign is for Therefor for the case of left-handed 
neutrinos we will have the dispersion relation, 

3 

E'^ =p^+ 2{r]i-r]2)^ (A. 6 ) 

Mpi 

and for antineutrinos we have, 

E^=p^-2{m-m)^ (A.7) 

Mpi 

We have dispersion relation for neutrinos and antineutrinos E"^ = p‘^ — 1 where 

^3 = — 2(71 — 72 ) for neutrinos and ^3 = 2(71 — 72 ) for antineutrinos. 

B Spinors Relation 

We assume that all the particles expect neutrinos follow the standard energy-momentum 
relation i.e, 

Ei = + (B.l) 
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where m* and pi are the mass and momentum of different particles (i = e,/r, r etc). The 
neutrinos follow the modified dispersion relation given in eq.(2.1). There exist very strin¬ 
gent bounds [22], which suggest that neutrino flavor is independent of their dispersion 
relation, so we assumed the universal dispersion relation for different flavor of neutrinos. 
We also define, 


^ ^ - 1 _ 

^ p ’ 


(B.2) 


where the function F{p) is the measure of the deviation of neutrino dispersion relation 
from the standard one [26]. In this framework, the modified Dirac equation for neutrino 
can be written as, 


{ij^do — iF{p)^ ■ d)'ip{x) = 0 


(B.3) 


where we have neglected the neutrino mass for simplification. Now we replace the Dirac 
field ijj in terms of the linear combination of plane waves i.e. 


V'(x) = u{p)e 


(B.4) 


using it, we get the following form of Dirac equation, 

{'y^E — F{p)^ ■ p)u{p) = 0. (B.5) 


Clearly, the positive energy solution of this equation will satisfy. 


E{p) = E{p)p, 


(B.6) 


we used these results in the derivation of the spinors sum of neutrinos, which comes. 


u%p)u%p) 

s=l,2 


I 0 p-a\ 
\p-a 0 j 


(B.7) 


where we assumed neutrino to be massless and defined p = {E,F{p)p). Following the Dirac 
algebra, we get the following result for spinor sum. 


Y u%p)u"{p) = = F{p)-i^^p^ (B.8) 

s=l,2 

where we used the result of eq.(B.6) for further simplification. For antiparticle when m = 0, 
there is an overall negative sign in eq.(B.5) and following the same procedure we obtain 
the same result, 

Y (B-9) 

s=l,2 
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